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Abstract 

Using the B s meson wave function extracted from non-leptonic B s decays, we evaluate the short 
distance contribution of rare decays B s — > l + £~ 7 (I = e, fi) in the standard model, including all 
the possible diagrams. We focus on the contribution from four-quark operators which are not taken 
into account properly in previous researches. We found that the contribution is large, leading to 
the branching ratio of B s — > being nearly enhanced by a factor 3 and up to 1.7 x 10~ 8 . The 

predictions for such processes can be tested in the LHC-b and B factories in near future. 

PACS numbers: 12.15.-y, 13.20.He 
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I. INTRODUCTION 



The standard model (SM) of electroweak interaction has been remarkably successful in 
describing physics below the Fermi scale and is in good agreement with the most experiment 
data. One of the most promising processes for probing the quark-flavor sector of the SM 
is the rare decays. These decays, induced by the flavor changing neutral currents (FCNC) 
which occur in the SM only at loop level, play an important role in the phenomenology of 
particle physics and in searching for the physics beyond the SM jl, 2|. The observation of 
the penguin-induced decay B — > X s j, B — > X s £ + £~(£ = e, /i) are in good agreement with 
the SM prediction, and the first evidence for the decay B s — > fi + fi~ was confirmed at the 
end of 2012 3|, putting strong constraints on its various extensions. Nevertheless, these 
processes are also important in determining the parameters of the SM and some hadronic 
parameters in QCD, such as the CKM matrix elements, the meson decay constant /b s , 
providing information on heavy meson wave functions. 

Thanks to the Large Hadron Collider (LHC) at CERN we have entered a new era of 
particle physics. In experimental side, in the current early phase of the LHC era, the 
exclusive modes such as B s — > £ + £~^ {£ = e, fi) are among the most promising decays due 
to their relative cleanliness and sensitivity to models beyond the SM jl|, Qj. In theoretically 
side, since no helicity suppression exists and large branching ratios as B s — > fi + n~ are 
expected. There are mainly two kinds of contributions of B s — > £ + £~^ in the SM: the short 
distance contribution which can be evaluated reliably by perturbation theory 4J while the 
ong distance QCD effects describing the neutral vector-meson resonances and J/\l/ family 
-[7] . As for the short distance contribution, it is thought in previous works that a necessary 
work is only attaching real photon to any charged internal and external lines in the Feynman 
diagrams of b — > s£ + £~ with statement that contributions from the attachment of photon 
to any charged internalpropagator are regraded as to be strongly suppressed and can be 
neglected safely , lSl, |9|, i.e., one can easily obtain the amplitude of B s — > £ + £ 7 by 
using the effective weak Hamiltonian of b — > s£ + £~ and the matrix elements ( , y\sOib\B s ) 
Oi = J^Pl, o-^ u q v PR directly. Therefore contributions from the attachment of real photon 
with magnetic-penguin vertex to any charged external lines are of course omitted l|, |2| or 



stated to be negligibly small 



9|. Another contribution from loop insertion of the lower order 



four-quark operators are also always neglected. We note that the complete contribution 
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seems to have been done in j^J, however it mainly concentrated on the long distance effects 
of the meson resonances, whereas the short distance contribution was indeed incompletely 
analyzed. A complete examination included all contribution to the processes in the SM is 
needed. 

As being well known, only short-distance contribution can be reliably predicted, and it is 
more important than the long-distance contribution from the resonances which is actually 
excluded partly by setting cuts in experimental measurements. Recently we showed that 
the contributions from the attachment of real photon with magnetic-penguin vertex to any 
charged external lines can enhance the branching ratios of B s — > by a factor about 2 

In this letter, we will extend our previous studies and use the B s meson wave function 
extracted from non-leptonic B s decays 11] to revaluate the short distance contribution from 



the all categories of diagrams of B s — > £ + £~7 decays. Special attention will be payed on the 
contribution from the four-quark operators, and a comparative study with previous work 
will be discussed. The paper are organized in the following, in sec. HI1 we analyse the full 
short distance contribution and present detail calculation of exclusive decays B s — > £ + £~^. 
The numerical results and the comparative study are given in sec. IHIt and the conclusions 
are given in sec. IIVI 



II. COMPLETE ANALYSIS ON SHORT DISTANCE CONTRIBUTIONS 



In order to simplify the decay amplitude for B s — > i + l 7, we have to utilize the B s 
meson wave function, which is not known from the first principal and model depended. 



Fortunately, many studies on non-leptonic B 
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13| and B s decays [ll| have constrained 



the wave function strictly. It was found that the wave function has form 
where the distribution amplitude <Pb s {x) can be expressed as fl^ |: 



(1) 



^>b s (x) = N Ba x 2 (l - x) 2 exp 



2 2 
m B x 



(2) 
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with x being the momentum fractions shared by s quark in B s meson. The normalization 
constant N Ba can be determined by comparing 

(0\&fy 5 b\B s ) = iN c f <p Bs {x)dxTr[Yl 5 {0 Bs +m Bs ) l5 }dx = -4N c tp Bg [ <f) Bs (x)dx (3) 
Jo Jo 

with N c being the number of quarks and 



(0\*fl B b\B a )=-if B .i% a 



(4) 



the B meson decay constant f Ba is thus determined by the condition 



<p Bs (x)dx = Jfi-fB.- 



(5) 




(a) 
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(d) 



FIG. 1: Feynman diagrams without effective vertex of b — > sj contribution to B s — > 77. 




FIG. 2: Feynman diagrams with effective vertex of b — > S7 contribution to B s — > 77 . The black 
dot stands for the magnetic-penguin operator O7 

Let us start with the quark level processes B s — > l + l~ r y which are subject to the QCD 
corrected effective weak Hamiltonian. The general effective Hamiltonian that describes 
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b — > s transition is given by 

AC 10 

74// = —jf v tbK E c jM°M> ( 6 ) 

where Oj (j = 1, . . . , 6) stands for the four-quark operators, and the forms and the corre- 



sponding Wilson coefficients Cj can be found in Ref. [15]. 

Generally, to describe all the short distance of the process B s — > £ + £~j, new effective 
operators for b — > S77 which are not included in ([6]) should be introduced. Corresponding 
feynman diagrams without and with effective vertex b — > 57 are shown in FIG. Q]and FIG. 
El respectively. When connect di-lepton line to one 7, operator b — > S77 may contribute to 
B s — > £ + £~ r y. Contributions from the such kind of diagrams with a photon attaching from 
internal charged lines to B s — > £ + £~ / -f are usually regraded as to be strongly suppresse d by 
a factor m^/m^ thus can be neglected safely 



2l, [8|, |9| . However, as pointed out in 



16] 



the conclusion of this is correct, but the explanation is not as what it is described. Here 
we address the reason more clearly: the contribution from diagrams FIG. [T] (a) and FIG. [TJ 
(b) are not suppressed. When applying effective vertex of b — > S7 to describe b —¥ S77 as 
shown in FIG. |2j internal quarks in the effective vertex are off-shell and such off-shell effects 
are also not suppressed. We have proven that the such two non-suppressed effects in FIG. 



[T]and FIG. [2] cancel each other exactly 16(. Therefore we can use the effective operators 
listed in Eq. ([6]) for on-shell quarks to calculate the total short distance contributions of 
B s -> £+£-7 in SM safely. 

The Feynmann diagrams contributing to B s — > £ + £~'-f at parton level can then be classified 
into three kinds as follows: 

1. Attaching a real photon to any charged external lines in the Feynman diagrams of 
b^s£ + £~; 

2. Attaching a virtual photon to any charged external lines in the Feynman diagrams of 
b — > S7 with virtual photon into lepton pairs; 

3. Attaching two photon to any charged external lines in the Feynman diagrams of four- 
quark operators with one of two photon into lepton pairs. 

Note that the third contribution is not considered in the previous studies except for Ref. 

which is the focus of this paper and will show the detail in the following. We also will 
discuss these contributions seperatly. 




(c) 



FIG. 3: Feynman diagrams that contribute to 
bution of O7, O9, O10 in tree level. 




(d) 



matrix elements B s — > £ + £ 7 with the contri- 



A. External real photon contributions 



The Feynman diagrams of the first kind of contributions are shown in FIG. [3j As the 
contribution from the FIG. |3] (c) and (d) with photon attached to external lepton lines, con- 
sidering the fact that (i) being a pseudoscalar meson, B s meson can only decay through axial 
current, so the magnetic penguin operator O7 's contribution vanishes; (ii) the contribution 
from operators O9, 0\q has the helicity suppression factor m^/m^, so for light lepton elec- 
tron and muon, we can neglect their contribution safely. These diagrams in FIG. |3] (a) and 
(b) are always regarded as the dominant contributions, and they have been considered by 



using the light cone sum rule jl, 2], the simple constituent quark model 



distribution amplitude extracted 
of B s — > £ + £~'~f at meson level as 

Ai = iN c Gee d 



rom non-leptonic B decays |9[. We rewrite the amplitude 



and the B meson 



m 



Pb s ■ k 

C w [C + ie a(3 ^p' Bs e l3 k u + C-PB s (e^k v - k^e v 



I5I 



(7) 
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The form factors in Eq. (J7J are found to be: 



9 ~ Z ^ °7 j 

r» - r e ffr n m b m Bs r eff r 

( - / 2 — ^9 < - y - — ^ 2 7 ~+ 



with the constant G = a em G fV^V^* / '(\/27r), and 



/ (i±-)«K(x)dx. 



The expression in Eq. ([7]) can be compared with Ref 



. a. 



(8) 



(9) 



^ — • JkAAAAAAA/ 
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(6) 




FIG. 4: Feynman diagrams of b — > sj with virtual photon into lepton pairs . 



B. External virtual photon contributions 



The Feynman diagrams of the second kind of contributions are shown in FIG. HI Con- 
tributions from the kind of diagrams are always neglected jl, 2] or stated to be negligibly 
small Note the B s meson wave functions used in this work and Ref. 9] are both from 
non-leptonic B s decays. However, as mentioned in the introduction the authors of Ref. {9] 
did not present the expression of the contribution from FIG. H] and only stated the con- 
tribution is numerical negligibly small. But such statement seems to be questionable, for 
that the pole of propagator of the charged line attached by photon may enhance the decay 
rate greatly which make some diagrams can not be neglected in the calculation. In these 
two diagrams, photon of the magnetic-penguin operator is real, thus its contribution to 
B s — > £ + £~ , y is different from the first kind contributions. We get the amplitude 10( : 

1 



A u = i2N c Gee d Cl 



,, ff m b m Bs 



T Pb s ■ q 



C+ [k,q ■ e — e^k • q — ie, ua ^h a q p ] [l-fl] , (10) 



with coefficients C + obtained by a replacement: 



C ~\- — C -\-{x — y x 

= N B I dx( 

J x - z - te 



x — z 

1 



y ->y = y 



ie) 



+ 



)x 2 (l 



X 



z — ie 



x) exp 



■in 



2u% 



B s 2 
X 



where z 



- and the first, second term in (fTTl) denotes the contribution from FIG. H] (a) 



and (b) respectively. Note that the contribution from FIG. [3] (a) is much larger than (b) 
since mg s <C (see next section) which is very easily understood in sample constituent 
quark model [8], i.e., 0b 3 (x) = S(x — m s /mB 3 ). However, the contributions from FigJH (a) 
and (b) are comparable, and pole in C + corresponds to the pole of the quark propagator 
when it is connected by the off-shell photon propagator. Thus the C + term may enhance 
the decay rate of B s — > £ + £~7 and its analytic expression reads 



2N Bs niz 2 ( y l - z) 2 exp 



+ N Be 
- N Bs 



2*4 



dx{ 

x + z 

i 



z 



1 + X 

1 



)x 2 (l + x) 2 exp 



1-i 

x 



— Z X 



dx 



ri- 



se 



exp 



in 



B s 2 
X 



2u 2 Bs 
2lj 2 Bs x 2 



(12) 




(a)0 3 ~ 6 





(6)0 3 ~ 6 
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(c)Oi~0 6 (q = c,b,s) 



S W\A<^ 

(d)Oi ~ 6 (g = c, b,s) 



FIG. 5: Feynman diagrams that contribute to the process B s — > £ + £ 7 with possible insertion of 
0\ to Oq in loop level . 
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C. Quark weak annihilation contributions 



Now we focus on the contributions from the diagrams for the four-quark operators which 
are not considered properly in previous works. The Feynman diagrams of the third kind of 
contribution are shown in FIG. The operator O7 is high order than four-quark operator 
0\ — 06, thus the contribution of loop diagram with operator 0\ — Oq insertion should be 
comparable with the tree level electro-weak penguin O7 — Oiq contributions listed above. To 
calculate the leading order matrix elements of b — > 37*7 shown in FIG. [51 we express the 
decay amplitude for b{p b ) —> s{p s )^*{ki)^{k 2 ) as 

A A Y 2 

A m (b -+ ^7*7) = i^V tb V; s ^s(p s )T» u b( Pb y(hy(k 2 ), (13) 

where pi, tS , k 12 denotes the momentum of quarks and photon respectively, e is the vector 
polarization of photon. We split the tensor T^ v into the momentum odd and even power 
terms for simplification. Keeping our physics goal in mind, Without loss of generality we 
assume that photon with momentum k\ is virtual and drop the terms proportional to k^ in 
the expressions. After straight calculation, we obtain 

T°J d (q) = e\ {-^ [ie va pxl X lsk«4 [k lfl f 2 (q) - k 2 J + (q)] + ie^ x l5 k^k lv f + {q) 

+ie fMVa ^ l5 [k?f + (q) - fc 2 a /-(g)]} , (14) 

+ *W75^/{(g)j(15) 



T%T{q) = -2^L { (k lu k 2 , - k x k 2 g, v 



im + (i - -)m 

z, 



"9 

where q denotes the quark in the internal line which two photons are attached and and e q 
is the number of electrical charge of the quark. The loop functions appearing in (|T5|) have 
forms as 

f±(q) = \ + — [ lng(z q ,Uq,x)— =F 7^- f lng(z q ,u q ,x), 

* Z q JO x ^ z q Jo 

h(q) = i|/ + (?) + /-(g)], / 2 (?) = ^-I/ + (g)-/-(g)], 
2u n r\ 



/s(g) = I ^g(z q ,u q ,x)dx, 

Z q Z a Jo 



H *q JO 

1 - (u + z)x(l - x) 

9{z,u,x) = — f r — , (16) 

1 — ux[ 1 — X) 



with /((g) = \ — /i(g), z q — 2k ^i 2 and u q — fy. Writing the amplitude in this ways, one 

)j (j = 1, . . . , 4) are inserted, only T° dc 

9 



can infer easily that for example, when operators Oj (j = 1, ... ,4) are inserted, only T° dd 



part can contribute to b — > 577 while the process receives both parts contributions as 0$$ 
are inserted. It is also easily obtained the similar result for the on-shell photons as in Ref. 



171 ] by setting u q = . 



With the amplitude of b — > 57*7 and B s wave function ready, we write the total contri- 
bution from FIG. 5 to exclusive decay of B s (pb s ) — > ^{k)i + i~ as 

A m = -2ie^J2 C ^ m ») [ T (PbM»K ~ + ^W^S.*^] [£<f£\, (17) 



<1 



with the form factors given by 



T{ = T 2 = T\ = Tl = 0; 
T\ = N C T% = N c e 2 fi(z c ,t); 



Kq=u,d,s,c,b J 



1 2 



T 5 



rp5 

1 2 



e 2 q h{z q ,t) + e 2 d N c [h(z b ,t)) + fr^t)} 

q=u,d,s,c,b 



N c I z b i/2 



4-Zk 1 

fs(z b ,t) + (l- T ^ t )(--f 1 (z b ,t)) 



(b s) 



■Nc E e 2 q f 1 (z q ,t)-2e 2 

q=u,d,s,c,b 



1 1 

- - fi(z b ,t)) + (b 5) 



Y e 2 q h{z q ,t)-2N c e 

q=u,d,s,c,b 



Z b i/a 2 
1 ,1 



z b i/2 2 



h{z h ,t)) + {b^s) 



(18) 



where q 2 in Eq. f TT7|) is the invariant mass square of lepton pair. The functions can be 
obtained directly from ( Fl6l) by redefined parameters z q = rn\ 2 /m Bg , t = q 2 /'m 2 Bs , 



A CM) 
f 3 (z,t) 



1 z f 1 dx r z — x(l — x) . 

2 + 1 - t 7 ~x~ ^z-tx(l -xy 

"dxln[-"- x{1 - x) 



f 



1 - t 



1 - t 



— te(l — sc) 



(19) 
(20) 



with explicit formula needed in calculation given by 



'dy 
2/ 



ln[l - vy(l - y) - ie] 



— 2arctan 2 a/ 43^ for u < 4; 

- 2Z7T In V^+v^=4 + 2 ln 2 v^+V^4 for ^ > 4 



/ dyln[l -vy(l - y) - ie] 
Jo 



11 — x 



1/2 



In 



iiT for x = 4/f < 1; 



2 arctan J— r for 2 = 4/w > 1. 



(21) 
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From Eq. (|T7|) it is clear that the contribution from four-quark operators to B s — > ^£ + £~ 
has the similar expression as that from magnetic-penguin operator with real photon to B s — > 
£ + £~-y. Thus the total matrix element for the decay B s — > £ + £~^ including contributions 
from three kinds of diagrams can be obtained easily by a shift to the form factors: 



ff _ 2m b m B . „ / x p^k\\ _ 2 Ps^kfs L y, 

q 2 V PB s -qJ\ q 2 e d J 

C 2 = C*"C_ - 2mb ^ C^C + ( 1 + + 2^-^ Y C 3 T{. (23) 

q 2 V q J q e d 

Finally, we get the differential decay width versus the photon energy J5 7 , 

- a - UF % b V t :\ 2 (m Bs - 2E,)E 1 [\C,\ 2 +\C 2 \ 2 + C 2 w (\C + \ 2 + |C_| 2 )] . (24) 



dE 7 108tt 



III. RESULTS AND DISCUSSIONS 



The decay branching ratios can be easily obtained by integrating over photon energy. In 
the numerical calculations, we use the following parameters 19| : 

a em = —, G F = 1.166 x l(T 5 GeV~ 2 , m b = 4.2GeV, 

\V tb \ = 0.999, \V U \ = 0.04, \V td \ = 0.0084 

m Bg = 5.37GeV, u Bs = 0.5, f Bs = 0.24GeV, r Bs = 1.47 x 10^ 12 s. 

m B o = 5.28GeV, cu Bd = 0.4, f Bd = 0.19GeV, r Bd = 1.53 x 10~ 12 s. 

The ratios of B s — > ^y£ + £~ are shown in Table [I] together with results of B^ s — > ^£ + £~ 
from this work and our previous research for comparison. The errors shown in the Table [J 
comes from the heavy meson wave function, by varying the parameter u Bi = 0.4 ±0.1, and 
lo Bs = 0.5 ± 0.1 Note that, the predicted branching ratios receive errors from many 
parameters, such as parameters meson decay constant, meson and quark masses. 

From the numerical results we conclude that unlike in decay B — > X .,77 where the four- 



quark operators just contribute a few percent to the branching ratio [17j, our numerical 
result shows that contribution from the four-quark operators to B s — > / ~f£ + £~ is large. It can 
be understood as follows: 

1. As pointed out in Ref. {9], the radiative leptonic decays are very sensitive probes in 
extracting the heavy meson wave functions; 

11 



TABLE I: Comparison of branching ratios in unit of 10 9 with previous calculations 



Branching Ratios 


Results 


(xl(T 9 ) 


The first kind of diaj 


grains 


The first two kind of diagrams 


Included all diagrams 


B s -»• j£+£~ 

B° d -> ji+e- 


5-i6l?:i 

0-2118& 


0.40^ 


1 7 Sfi +4 - 55 
u.oo_ 12 



2. Values of the Wilson coefficients Cj(mb) (j = 3, . . . , 6) are at order of 10~ 2 , indicating 
that contribution to the corresponding operators via T J is less important compared 
with those from Oi )2 ; 

3. From Eq. ( II 8p . one can infer easily that the four-quarks contribution to the form fac- 
tors in ( 123]) have coefficient (N c Ci+C2)T 2 fsjed while the contribution from magnetic- 
penguin operator with real photon, mbmsj \pb s • q)C^ • C+- Note {N c Ci + C2) / e<i and 
C7 can be comparable and have the same sign in C\ and opposite sign in C 2 . How- 
ever, with T( = for j = 1,2 thus comparable contribution studied in this work and 



in Ref. 



10( is expected, leading to enhancement of branching ratios of B s — > r )£ + £ 



when new diagrams are taken into account. 

4. The predicted short- distance contributions from quark weak annihilation as well as 
the magnetic-penguin operator with real photon to the exclusive decay are large, and 
the branching ratios of B s — >■ l + l^ r y are enhanced nearly by a factor 3 compared with 
that only contribution from magnetic-penguin operator with virtual photon and up to 
1.7 x 1CT 8 , implying the search of B s — > £ + £~ / ~y can be achieved in near future. 

5. Due to the large contributions from magnetic-penguin operator with real photon and 
quark weak annihilation , the form factors for matrix elements (7157^(1 — , ~f^)b\B s ) 
and ^scr^l ±^ 5 )q u b\B s ) as a function of dilepton mass squared q 2 are complex and 



not as simple as l/(q 2 — q^) 2 where g 2 is constant 18] . The B s — > 7 transition form 



-7|. For 
f tensor 

and pseudotensor currents with direct emission of the virtual photon from quarks are 



factors predicted in this works have also some differences from those in Ref. 
instance, Ref. \&\ predicted the form factors F TV (q 2 ,0), F T A(q 2 ,0) induced by tensor 
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only equal at maximum photon energy, whereas the corresponding formula in this 
work have the same expression as — ed p^ m ^ s C+ oc l/(q |2 — Qo) in Eq. (J7J). Furthermore, 
the form factors are larger than previous predictions. 

To clarify things more clear, we think it is necessary to present a few more comments 
about the calculation of Ref. [5|, as mentioned in introduction. In order to estimate the 
contribution of direct emission of the real photon from quarks, the authors of Ref. 5| calcu- 
lated the form factors Fta,tv{®, Q 2 ) by including the short- distance contribution in q 2 — > 
limit and additional long-distance contribution from the resonances of vector mesons such as 
p°, u° for Bd decay and <fi for B s decay. Obviously, this means the short- distance contribu- 
tions were not appropriately taken into account. Moreover, if Fta,tv(®, <f) = ^ta,tv(0, 0) 
stands for the short distance contribution, it seems to double counting since in this case 
photons emitted from magnetic-penguin vertex and quark lines directly are not able to be 
distinguished. 

We also note that for contribution from the weak annihilation the authors of Ref. 
only took into account u and c quarks in the loop by axial anomaly as the long distance 
contribution, they concluded that the anomalous contribution is suppressed by a power of 
a heavy quark mass. We believe that only anomalous contribution to account contribution 
from the weak annihilation is insufficient. Our numerical result shows that the contributions 
from weak-annihilation diagrams are large and can not be neglected. 



IV. CONCLUSION 

In summary, we evaluated short distance calculation of the rare decays B s — > r yl + l~ in the 
SM, including contributions from all four kinds of diagrams. We focus on the contribution 
from four-quark operators which are not taken into account properly in previous researches. 
We found that the contributions are large, leading to the branching ratio of B s — > £ + £~7 
being nearly enhanced by a factor 3. In the current early phase of the LHC era, the exclusive 
modes with muon final states are among the most promising decays. Although there are 
some theoretical challenges in calculation of the hadronic form factors and non-factorable 
corrections, with the predicted branching ratio at order of 10~ 8 , B s — > fi + fi~ / -f can be 
expected as the next goal after B s — > fi + fi~ since the final states can be identified easily 
and the branching ratios are large. Experimentally, B s — > / u + /i~7 mode is one of the main 

13 



backgrounds to B s — > and thus it is already taken into account in B s — > ji + n~ 

searches jjj. Our predictions for such processes can be tested in the LHC-b and B factories 
in near future. 

Acknowledgment s 

This work was supported in part by the NSFC No. 11005006, 11172008. 



[1] T.M. Aliev, A. Ozpineci, M. Savci, Phys. Rev. D 55,7059 (1997). 

[2] Z. Heng, R. J. Oakes, W. Wang, Z.H. Xiong and J. M. Yang, Phys. Rev. D 77, 095012 
(2008); G. Lu, Z.H. Xiong, Y. Cao, Nucl. Phys. B 487, 43 (1997); Z.H. Xiong, J. M. Yang, 
Nucl. Phys. B 602, 289 (2001) and Nucl.Phys. B 628, 193 (2002); G. Lu et. al. , Phys. Rev. 
D 54, 5647 (1996); Z. H. Xiong, High Energy Phys. Nucl. Phys. 30, 284 (2006). 

[3] RAaij et al. [LHCb Collaboration], Phys. Rev. Lett. 110, 021801 (2013) [ arXiv:1211.2674l 
[Unknown]]. 

[4] G. Buchalla, A. J. Buras, M. E. Lautenbacher, Rev. Mod. Phys. Vol 68, No. 4, 1125. 
[5] D. Melikhov and N. Nikitin, Phys. Rev. D70, 114028, (2004). 
[6] F.Kruger and D.Melikhov, Phys.Rev. D67, 034002 (2003); 

[7] N. Nikitin, I. Balakireva, and D. Melikhov. llarXiv: 1101 .4276^ 1 [hep-ph]; C.Q. Geng, C.C. Lih, 
W.M. Zhang, Phys.Rev. D62, 074017 (2000); Y. Dincer and L.M. Sehgal, Phys.Lett. B521, 
7 (2001); S.Descotes-Genon, C.T. Sachrajda, Phys.Lett. B557, 213 (2003). 

[8] G. Eilam, C.-D. Lii and D.-X. Zhang, Phys. Lett. B 391, 461 (1997); G. Erkol and G. Turan, 
Phys. Rev. D 65, 094029 (2002). 

[9] J.-X. Chen, ZY. Hou , C.-D. Lii, Commun. Theor. Phys. 47, 299 (2007). 
[10] W. Wang, Z. -H. Xiong and S. -H. Zhou. l arXiv:1207.1978l [hep-ph], 

[11] Y. Li, CD. Lii, Z.J. Xiao, and X.Q. Yu, Phys. Rev. D 70, 034009 (2004); X.Q Yu, Y. Li and 
CD. Lii, Phys. Rev. D 71, 074026 (2005) and 2006 Phys. Rev. D 73, 017501. 

[12] C.-D. Lii, K. Ukai and M.-Z. Yang, Phys. Rev. D 63, 074009 (2001); Y.-Y. Keum, H.-n. Li 
and A.I. Sanda, Phys. Lett. B 504, 6 (2001), Phys. Rev. D 63, 054008 (2001); H.-n. Li, 2001 
Phys. Rev. D 64, 014019 (2001); S. Mishima, Phys. Lett. B 521, 252 (2001); C.-H. Chen, 



11 



Y.-Y. Keum, and H.-n. Li, Phys. Rev. D64, 112002 (2001) and Phys. Rev. D 66, 054013 
(2002). 

[13] CD. Lii, Eur. Phys. J. C 24, 121 (2002) and Phys. Rev. D 68, 097502 (2003); Y.-Y. Keum 
and A.I. Sanda, Phys. Rev. D 67, 054009 (2003); Y.-Y. Keum, etaL, Phys. Rev. D 69, 094018 
(2004); J. Zhu, Y.L. Shen and CD. Lii, 2005 Phys. Rev. D 72, 054015, and Eur. Phys. J. 
C41, 311 (2005); Y. Li and CD. Lii, Phys. Rev. D73, 014024 (2006); CD. Lii, M. Matsumori, 
A.I. Sanda, and M.Z. Yang, Phys. Rev. D 72, 094005 (2005). 

[14] C.-D. Lii and M.-Z. Yang, 2003 Eur. Phys. J. C 28, 515. 

[15] M. Misiak, Nucl. Phys. B 393, 23 (1993). 

[16] F.L. Dong, WY. Wang, and Z.H. Xiong, Chin. Phys. C 35, 6 (2011). 

[17] J.J. Cao, Z.J. Xiao and G.R. Lu, Phys. Rev. D 64, 014012 (2001); L. Reina, G. Ricciardi, 

and A. Soni, Phys. Rev. D 56, 5805 (1997). 
[18] G. Eilam, I. Halperin and R. R. Mendel, Phys. Lett. B 361, 137 (1995). 
[19] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010) and 2011 partial 

update for the 2012 edition. 



15 



